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CONSEQUENCES OF COSMIC MICROWAVE BACKGROUND-REGULATED STAR FORMATION 
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ABSTRACT 

It has been hypothesized that the cosmic microwave background (CMB) provides a temperature 
floor for collapsing protostars that can regulate the process of star formation and result in a top-heavy 
initial mass function at high metallicity and high redshift. We examine whether this hypothesis has 
any testable observational consequences. First we determine, using a set of hydrodynamic galaxy 
formation simulations, that the CMB temperature floor would have influenced the majority of stars 
formed at redshifts between z = 3 and 6, and probably even to higher redshift. Five signatures 
of CMB-regulated star formation are: (1) a higher supernova rate than currently predicted at high 
redshift; (2) a systematic discrepancy between direct and indirect measurements of the high redshift 
star formation rate; (3) a lack of surviving globular clusters that formed at high metallicity and 
high redshift; (4) a more rapid rise in the metallicity of cosmic gas than is predicted by current 
simulations; and (5) an enhancement in the abundances of a elements such as O and Mg at metallicities 
—2 < [Fe/H] < —0.5. Observations are not presently able to cither confirm or rule out the presence 
of these signatures. However, if correct, the top-heavy IMF of high-redshift high-metallicity globular 
clusters could provide an explanation for the observed bimodality of their metallicity distribution. 
Subject headings: stars: formation — cosmic microwave background — stars: luminosity function, 
mass function — galaxies: stellar content — globular clusters: general — stars: 
abundances 



1. INTRODUCTION 

The initial mass function (IMF) describes the relative 
numbers of stars formed with different masses. Observa- 
tions suggest that the IMF has a power-law form at the 
high mass end {M > 1 Mq , with a slope near the c anon- 
ical Salpeter value of a = —2. 35: lSalpeteiill955D. and 
a turnover at the low-mass end (IKrouDall200inChabrieri 
[2nf)3h . 

The detailed form of the IMF is important because 
stars of different mass have different mass-to-light ratios, 
different lifetimes, and different effects on their surround- 
ings. Therefore, the luminosity, chemical enrichment, 
and energetic feedback due to a stellar population, as 
well as the evolution of these quantities with time, all 
depend on the IMF. 

In the local universe, observations indicate that the 
IMF is universal, showing no eviden ce for variatio n be- 
tween different star formation events ()Kroupall200lD . Al- 
though we h ave no full theory that ex plains the origin of 
the IMF fsee lMcKee fc Ostrikeii[2007l for a good review) , 
it has recently been discovered that the functional form 
of the IMF is the same as for the mass function of preste l- 
lar cores, but shifted to lower mass (|Alves et aLlIM)! . 
This suggests that each star forms with a constant frac- 
tion of its core mass, and the functional form of the IMF 
is set by the process of gas fragmenting into cores. As 
such, it should be related to how the Jeans mass evolves 
within a collapsing protostellar environment. This leads 
to the intriguing possibility that star formation in en- 
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vironments where the thermal and density evolution is 
dramatically different than in the local universe could 
result in different IMFs. 

One such environment is the virtually metal-free gas 
that formed the very first "Population III" stars. With- 
out any metals, cooling below 10"* K is very inefficient, 
and the fragmentation mass of the gas remains high. De- 
tailed hydrodynamic simulations of the formation of Pop- 
ulation III stars suggest that stars formed from such zero- 
metallicity gas follow a very top-heavy IMF, with a me- 
dian stellar mass perhaps as high as 100 Mq (|Abel et al.l 
I200I . A more normal IMF is obtained when the gas 
reaches a criti cal value ZrT\t , somewhere between 10~^ 
and 10-3 Zq (|Bromm et all 120011: iSmith fc Sigurdssonl 
l20?il . 

Another important but less well-studied regime in 
which a top-heavy IMF has been proposed is for high- 
redshift /ii!(//i-metallicity gas. Unlike primordial gas, 
which cannot cool efficiently enough to fragment, high 
metallicity gas cools very efficiently. However, at high 
redshift the cooling (and therefore fragmentation) comes 
to an abrupt halt when it reaches the temperature of 
the cosmic microwave background radiation (CMB). Be- 
cause it is thermodynamically impossible for the proto- 
star to cool below the CMB via radiative mechanisms, 
and cooling via adiabatic expansion is not relevant to a 
collapsing protostar, the CMB sets a temperature floor. 
This abrupt halt to fragmentation results in a top-heavy 
IM F. This argum ent was put forward in analytic form 
bv lLarsonI (I2OO50. and was fu rther followed up numer- 
icallv bv lOmukai et al.l (|2OO50 . who used a single-zone 
collapsing protostell ar model wi t h an advanced chemi- 
cal network, and by ISmith et al.l ([20091 ) , who performed 
ab initio 3D hydrodynamic simulations of high redshift 
star formation. These latter authors, in particular, argue 
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that high-redshift high-metallicity star formation must 
produce very massive stars due to the influence of the 
CMB. 

The idea that the CMB regulates star formation is in- 
triguing and could be very important. However, little 
work has yet been done on estimating whether it would 
have had any influence on real star formation, or on de- 
termining the observational consequences of such regula- 
tion. The goals of the present paper are to estimate the 
regimes in which this effect could be important, and to 
examine possible consequences that could be confirmed 
or ruled out by observations of the local universe and at 
high redshift. In § [21 we provide a relationship between 
the redshift of star formation and the critical metallicity 
for CMB regulation. In §[3l we use high resolution cosmo- 
logical hydrodynamic simulations to estimate the amount 
of star formation that occurred at different metallicities 
as a function of redshift. We bring these together in § |4] 
and find that, if the CMB does couple to star-forming 
gas, the majority of high-redshift star formation must 
have happened in the CMB-regulated regime, and exam- 
ine the effects of this regulation on supernova and gamma 
ray burst rates, direct vs. indirect measurements of the 
high redshift star formation rate, the metallicity distri- 
bution of old globular clusters, and the cosmic evolution 
of both global metallicity and a abundances. Finally, we 
present our conclusions in § [5] 

2. ESTIMATING Zcmb 

According to the CMB-regulation hypothesis, star- 
forming gas that is sufficiently metal-rich cools rapidly 
to the CMB temperature, whi ch halts fragmentation and 
resul ts in a top-heavy IMF (jLarso ^ 120051 : ISmith et al.l 
l2009f) . The metallicity is a key parameter because the 
cooling rate at T < lO'* K is completely dominated by 
meta l line emission for all but the lowest metallicities 
(e.g. ISmith et alll2008D . There should therefore be at 
each redshift z a critical metallicity, Zqmb,(z\ above 
which CMB regulation is important. 

The best estim ates of .^cmb to date come from 
ISmith et al . (2009), who performed ah initio simulations 
with three different sets of initial conditions and a vari- 
ety of metallicities. The amount of fragmentation they 
found in the simulations depended on metallicity: at low 
metallicities, no fragmentation was seen because the gas 
remained hot; at intermediate metallicities, the gas frag- 
mented; and at high metallicities, there was again no 
fragmentation because the gas cooled so quickly that it 
reached the CMB temperature. Because simulations can 
only be performed for a discrete set of metallicities, we 
cannot determine the exact value of Zqm^ in each case, 
but we can bound it on the lower end by the most metal- 
rich simulation where fragmentation occured, and on the 
upper end up the least metal-rich simulation where frag- 
mentation did not occur. These are plotted in Figure[T]as 
the connected diamonds, at the redshifts at which each 
simulation collapsed. Note that although the Set 2 sim- 
ulations collapse at a redshift intermediate between the 
Set 1 and Set 3 simulations, their implied Zqmb is higher, 
demonstrating that the boundary between the normal 
and CMB-regulated regimes is not purely a function of 
metallicity, but also depends on t he details of the initial 
conditions; 'Japps en et al.l ()2009D make a similar point 
regarding the lower metallicity threshold, ZcHt- 
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Figure 1. Estimated evolution of Zqmb as a function of red- 
shift (solid line), above which fragmentation is inhibited due to 
the CMB. The dot-dashed line shows the critical metallicity Z^^n, 
below which fragmentation is inhibited due to inefficient cooling; 
we use the v alue log Z/Zq = —3.9, from the Set 2 simulations of 
ISmith et aLl 12009). The region in between these lines, where a nor- 
mal IMF is expected, is shaded. The dotted vertical line at 2 = 2.7 
denotes when the CMB temperature drops below 10 K; the CMB 
is not expec ted to be important below this point. The results of 
ISmith et aLl (2009) are also shown; the squares denote the values 
of Zcrit they calculated for each of their sets of initial conditions, 
while the pairs of connected diamonds denote the range of possi- 
ble values of Zqmb ■ Note that because the same initial conditions 
collapsed at different redshifts for different metallicities, the pairs 
of diamonds are not vertical. 

It is difficult to determine the consequences of CMB 
regulation purely from these simulations because very 
little star formation occurs at such high redshifts. We 
must therefore devise a method of estimating the evo- 
lution of Zqmb as a function of redshift . To do this, 
we no te th at in the numerica l models of lOmukai et al] 
([2005h and ISmith et~all (pOOl . the CMB effectively acts 
as a temperature floor to the cooling collapsing proto- 
stellar core. Without the CMB, the normal evolution 
of the core is to cool and collapse until heat is released 
by rapid H2 formation via three-body collisions, which 
drives a rise in temperature. Once all of the hydrogen 
is in molecular form, cooling once again dominates and 
the temperature drops. Eventually the core becomes op- 
tically thick, at which point the dust emission can no 
longer cool the core, and the temperature rises again due 
to compressional heating. There is a minimum temper- 
ature Tniin(Z) that the gas reaches during this evolution 
(either immediately before the onset of the three-body re- 
action or immediately before the core becomes optically 
thick), which depends on its cooling history and there- 
fore its metallicity. Our ansatz is as follows: at redshift 
z, the CMB regulates star formation in gas of metallicity 
Z if the CMB temperature 

7cMb(^) > T^in{Z), (1) 

where the CMB temperature at redshift z is 

Tcmb{z)^To{1 + z), (2) 

for present-day CMB temperature Tq = 2.726 K. If 
the effect of the CMB is to act as a temperature floor, 
this condition defines the regime where the CMB can 
ha ve an effect. We det ermine r,„in(Z) from the tracks 
of lOmukai et al.l ()2005f ). who modelled the evolution 
of temperature and density for a series of metallicities 
Z/Zq = 0, 10-^ 10-5, 10-^ 10-3, 10-2, 10"i, 10°. For 
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Figure 2. Density-temperature evolution of a protostellar core 
witii rnetallicit y Z = 10~^ Zq (10~^ Zq) in red (blue) from 
[Omukai ct al] ll|2005l ). with a temperature floor set by the CMB 
temperature at each of the four labelled redshifts. The horizon- 
tal dashed lines denote the CMB temperature at each redshift. 
For Z = 10-2 Zq, Tjnin = 14.2 K, while for Z = 10"^ Zq, 
Tlnin = 34.1 K. The diagonal dashed gray lines indicate lines of 
constant Jeans mass, and are labelled in Mq. 

each of these metalhcities, we can then determine the 
redshift where equation ([l} is satisfied. 

This is demonstrated in Figure [2l where we have 
plotted the density-temperature evolution of protostellar 
cores of two differe i it met alhcities {log Z/Zq = —2, —3) 
from lOmukai et aD (|2005f ). In each of the four panels, we 
have imposed a temperature floor given by Tcmb at a dif- 
ferent redshift. At z = 3, Tcmb — 10.9 K, which is colder 
than either core ever reaches, and the CMB has no effect. 
At z = 4.2, Tcmb = 14.2 K, which is the exact minimum 
temperature that the \og Z/Zq = —2 protostar reaches; 
if the CMB temperature were any higher, it would af- 
fect the core evolution. This can be seen at z = 11.5, 
when Tcmb = 34.1 K; the logZ/Z© = —2 core does 
not cool nearly as far, and the CMB temperature just 
reaches the minimum temperature of the log Z/ Zq = —3 
protostar. At z = 15, the temperature floor imposed 
by the CMB clearly affects the evolution of both cores. 
We therefore adopt ZcmbIz = 4.2) = 10^^ Zq and 
Zcmb{z = 11.5) = 10"^ Zq. We pe rform the same calcu - 
lation for each metallicity track in lOmukai et al.l ()2005f ) 
and plot the relationship as the solid line in Figure [T] 

Fragmentation is only expected to be efficient when the 
Jeans mass rapidly decreases, i.e. when the tracks move 
downward. Therefore the characteristic mass scale is the 
Jeans mass when the core stops cooling (Larson 2005]). 
Lines of constant Jeans mass are shown as the diagonal 
dashed gray lines in Figure [2J For the Z = 10 Zq 
track, this characteristic scale is 0.2 Mq, but at redshift 
z = 11.5 it rises to 50 Mq. The top-heavy IMF is a direct 
consequence of this dramatic change in the characteristic 
mass scale at the end of fragmentation. 

Because the temperature of the protostellar core only 
just reaches the CMB temperature, it is unlikely that 
the CMB has much effect at the exact redshift we calcu- 
late. The Tcmb(^), and therefore z, that we calculate for 
each metallicity must be a slight underestimate, or corre- 
spondingly the ^cMB that we calculate at each redshift 
must be a slight underestimate. 

Our estimated Zcm b matches the simulation results 
of iSmith et"all ()2009[ ) very weU for their Set 1 initial 
conditions, mildly underpredicts it for their Set 3 intial 
conditions, and underpredicts it by about an order of 



magnitude for their Set 2 initial conditions. We should 
therefore expect that the transition between a normal 
and top-heavy IMF occurs at metalhcities somewhere be- 
tween our predicted Zcmb and a metallicity ten times 
larger, with the details depending on the properties and 
formation history of the individual halos in which the 
star formation occurs. 

Although our estimate for Zqub is well-defined down 
to z « 0.2, it is unlikely that CMB regula tion is truly 
important at these redshifts. As noted by ISmith et al.l 
(j2009), gas clouds in the local ISM are not observed to 
cool below 10 K; therefore, when the CMB drops be- 
low this temperature, it can no longer affect the thermal 
evolution of protostellar gas. This occurs at z = 2.7, 
denoted by the vertical dotted line in Figure [TJ 

As discussed earlier, a top-heavy IMF is also expected 
for Population III stars, with metalhcities below Za-it, 
denoted by the horizontal dot-dashed line in Figure [TJ A 
normal IMF is therefore expected in the shaded region 
of Figure [l] at lower metalhcities, cooling is too ineffi- 
cient for fragmentation to occur, while at higher metal- 
hcities, the gas cools immediately to the CMB tempera- 
ture where it becomes isothermal and does not fragment. 
Although the estimates of the boundaries of the region 
are rough, they provide a guideline for the star formation 
events that would have been influenced by the CMB. 

It is interesting that our derived Zqmb rises substan- 
tially from < 10-3 Zq at z > 10 to ~ lO-^-^ at z ~ 3. 
Such metalhcities are typical for old stellar populations, 
and it is therefore plausible that CMB regulation could 
have been important for stars formed at these redshifts. 

3. COSMIC METALLICITY EVOLUTION 

In order to determine the effects of CMB-regulated star 
formation, we must estimate, at each redshift z, the frac- 
tion of stars that formed with metalhcities Z > Zcmb(-z)- 
Observational measurements are available at z < 3 from 
stellar population modelling of the integrated spectra 
of galaxies from t he Sloan Digital Sky Survey (SDSS; 
iPanter et al.ll2008f ). but, as discussed in § [21 the effects 
of the CMB are only likely to be significant at z > 3. 
We must therefore use theoretical models to estimate the 
cosmic evolution of the metallicity of star-forming gas. 

3.1. Simulations 

Our simulations were performed as part of the Mc- 
Master Unbiased Galaxy Simulations project (MUGS), 
a campaign to construct high resolution simulations of a 
large set of T* galaxies that randomly sample the sites of 
galaxy formation, including a full range of modern galaxy 
formation physics. Full details of the MUGS simulations 
will be presented in Stinson et al. (in preparation); we 
provide an overview of the most important properties of 
the simulations below. 

The simulations were performed using a WMAP 3 
ACDM cosmology with Hq = 73 km s-^ Mpc-\ Qm = 
0.24, n^ = 0-76^ f^bary = 0.04, and erg = 0.76 
([Spergel et al.ll2007f ). Halos were chosen from a uniform- 
resolution 256'^ dark matter-only simulation in a box of 
side 50 fe-^ Mpc using the friends-of-friends algorithm 
(jPavis et al.l 119851 ). A random selection of isolated ha- 
los with masses 4 x 10" M© < M < 2 x 10^^ 
were chosen for resimulation at higher resolution with 
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full baryonic physics. The highest resolution dark mat- 
ter region, encompassing all matter out to 5 Tvir, has 
particle mass mDM = 1-1065 x 10^ Mq, while gas par- 
ticles are placed inside 3 r^ir and have initial masses 
mgas = 2.2131 X 10^ Mq. There are typically ~ 1000000 
gas and high resolution dark matter particles each in the 
refined regions of the resimulations, with the exact num- 
ber depending on the halo mass and the geometry of the 
Lagrangian region that collapses into the z = halo. 

The simulations were evolved using the parallel 
SPH code GASOLINE (Wadsley ct al. 2004). gaso- 
line solves the equations of hydrodynamics using SPH 
and self-gravity usin g the Barnes-Hut tree algorithm 
(|Barnes fc Hutl[l986l ). and includes radiative cooling, an 
ultraviolet (UV) background, star formation, and ener- 
getic and chemical feedback. 

The cooling is calculated from the contributions of 
both primordial gas and metals as Atotiz, p,T, Z) = 
AHi,Hoi,Hcii(^;,/3,r) + -^A^ctiii.Zeiz,p,T). The first 

term employs atomic cooling based on a gas with pri- 
mordial composition heated by a uniform UV ionizing 
backgroun d, adopted from Haardt & Madau (in prepa- 
ration; see lHaardt fc Madaull200lD with ra t es coe fficient 
closely matching those cited in lAbel et al.l (|1997| ) , while 
the metal cooling grid is constructe d using Cloudy (ver- 
sion 07.02, last described by Ferla nd et al.l ( 1998D'). as- 
sumin g ionization equilibrium, as described in Shen et aLl 
()2009( ). The UV background is used in order to calcu- 
late the metal cooling rates self-consistently. The cooling 
lookup table is linearly interpolated in three dimensions 
(i.e., p, z, T) and scaled linearly with metallicity. 

The star formation and feedback recipes are based 
on the "bla s twave model" described in detail in 
IStinson et al.l (|2006f) . but with the addition of clus- 
tered supernovae to account for the clustered nature 
of star formation. Star formation can occur in gas 
particles that are dense (rimin = 0.1 cm~^) a nd cool 
(Tmax = 15,000 K), calibrated to match the Ke nnicutO 
(ll998D Schmidt Law fo r the Isolated Model Milky Way 
in IStinson eFall ()2006D . 

At the resolution of these simulations, each star parti- 
cle represents a large number of stars (6.32 x 10'* Mq). 
Thus, each particle has its stars partitioned into mass 
bins based o n the initial mass function presented in 
iKroupa et al.l (|1993D . These masses are correlat ed to 
stellar lifetimes as described in lRaiteri et al.l ()1996D . We 
stochastically determine when a star particle releases 
feedback energy so that a minimum of 30 supernovae 
worth of energy is released concurrently to reflect the 
clustered nature of star formation. The explosion of these 
stars is treat ed using the analyt i c mod el for blastwaves 
presente d in IMcKee fc Ostrikeri (jl977D as described in 
detail in IStinson et al.l (j2006D . While the blast radius 
is calculated using the full energy output of the super- 
nova, less than half of that energy is transferred to the 
surrounding ISM, Esn = 4 x 10^° ergs. The rest of 
the supernova energy is radiated away. Iron and oxygen 
ar e produced in SNII a ccording to the analytic fits used 
in iRaiteri et al.l (|1996( ). The iron and oxygen are dis- 
tributed to the same gas within the blast radius as is the 
supernova energy ejected from SNII. E ach SNIa produces 
0.63 MfT) iron and 0.13 Mq oxygen ()Thielemann et aTl 
I1986D and it is ejected into the nearest gas particle for 
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Figure 3. Metallicity of stars as a function of their forma- 
tion redshift. The solid line is the mean from the MUGS sim- 
ulations, the dotted line is the mean from the simulation of 
Tpave & QppfS ( 2007), and the dot-dashed line is the ob- 

served mean metallicity of star-forming gas at each redshift, as 
determined from spectral synthesis modelling of SDSS galaxies 
CPanter et al.,.200ai) . 

SNIa. 

We have implemented diffusion of all scalar SPH quan- 
tities, particularly metal content and thermal energy, as 
described in Shen et al. (2009), which is required to cor- 
rectly model even simple process es such as convection 
and Rayleigh- Taylor instabilities (jWadslev et al.l I2008[) 
and to account for mixing in turbulent outflows. 

MUGS galaxies are labelled by their group number 
in the list returned by the friends-of-friends algorithm. 
The simulations analyzed in this work are MUGS gl536, 
g5664, g7124, gl5784, g21647, g22437, g22795, and 
g24334. Only stars within the virial radius of the main 
galaxy are considered. 

3.2. Metallicity Evolution 

The mean metallicity of stars formed in the simulation 
is shown in Figure [3] as a function of their formation red- 
shift. The solid line shows the results for all stars within 
MUGS simulated galaxies, and shows that stellar nietal- 
licities rise from ~ —3, for those formed at z > 10, to 
nearly solar for those formed aX z = Q. Our simulations 
are of L* galaxies, while the majority of star formation 
at high redshift occurred in larger galaxies, which formed 
their metals earlier than less massive galaxies. Our deter- 
minations may therefore underestimate the metallicities 
of a universal sample of stars formed at high redshift. 

We have confirmed that our results are consistent with 
those obtained with completely different codes: the dot- 
ted line, which shows the star-formation-rate-weighted 
mean metallicity of gas from the GADGET2 simulations 
of Dave & Oppenhcimer (2007) (the "SFR-wcighted" 
line in their figure 2) and should be directly compara- 
ble, shows reasonable agreement with our results over 
the entire range < z < 6 that they plotted. Where 
our re sults deviate from those of iDave fc Oppenheimeii 
(|200l . it is in the sense that the MUGS metallicities are 
lower. 

The best observational measurements of stellar metal- 
licity as a f u nction of formation redshift come from 
iPanter et al.l (|2008l ). who performed spectral synthesis 
modelling of SDSS galaxies at redshifts ranging from 0.1 
to 3. The mean metallicities of stars inferred to have 
formed at each redshift is shown as the dot-dashed line 
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Figure 4. (Top) The data points denote tlie median metallicity 
of stars formed at eacii redshift within the MUGS simulations, 
while the error bars denote the 10th and 90th percentiles. Stars 
formed within the dark shaded region, which is identical to the 
shaded region in Figure [T] are expected to form with a normal 
IMF. The light shaded region extends to metallicities 1 dcx higher 
than Zqmb, where stars may form with either a normal or top- 
heavy IMF, if CMB regulation occurs. All stars formed in the 
white region are expected to exhibit a top-heavy IMF if the CMB 
regulation hypothesis is correct. (Bottom) Fraction of stars formed 
with metallicities greater than ZqmBj i-e. within the light shaded 
or white regions of the top panel. The dashed line denotes the 
fraction of stars formed with metallicities greater than 10 Zqmb, 
i.e. within the white region of the top panel. 

in Figure [31 Unlike the simulation predictions, the ob- 
servations show essentially no drop in stellar metallicity 
out to z = 3. This is mainly because the total stel- 
lar mass is dominated by the most massive galaxies, 
which formed most of their stars very early; however, 
the stellar metallicities within galaxies of stellar mass 
3 X 10^° Mq < < 1 X 10" Mq show almost identical 
behavior, and this is very nearly the same mass range as 
the MUGS galaxies {4x10^° Mq < M, < 1.2xlO"M0), 
and the mass range that includes the Milky Way (M* 
5 X 101° Mq; iKlvpin et al.l[200l IWidrow et an[2051 . 
We may therefore conclude that our simulations do not 
overpredict the metallicity of star-forming gas at high 
redshifts, and may even underpredict it, an issue we will 
return to in § 14.51 

4. CONSEQUENCES OF CMB REGULATION 
4.1. When was the CMB important? 

We have plotted our estimate of the metallicity of star- 
forming gas from § [3] over our estimate of the critical 
metallicity for CMB regulation, Zcmb, from § [H in the 
top panel of Figure [H Note that unlike in Figure |31 here 
we plot the median metallicity, along with the 10th and 
90th percentiles. 

The first remarkable thing to note about Figure |4] is 
that the median star was formed with Z > Zqmb at al- 
most all redshifts greater than 2.7. In other words, if the 
CMB regulation hypothesis is correct, the majority of 
stars at high redshift were formed in the CMB-regulated 
regime. While this statement is subject to uncertainties 
in both our estimate of Zcmb, and in the cosmic metal- 
licity evolution measured from the simulations, it would 
appear that the CMB influenced a significant amount of 
star formation at high redshift. 

In the bottom panel of Figure IH we have plotted the 
fraction of stars that form aX Z > Zqmb in the sim- 



ulation. As expected from the top panel, this is very 
high: over 80% at virtually all redshifts where the CMB 
temperature is higher than 10 K. Given the uncertain- 
ties in both the estimate of Zcmb and in the metallicity 
evolution predictions of the simulation, we have also de- 
termined the fraction of stars formed with metallicity 
Z > 10 Zcmb (dashed line). We consider it unlikely that 
the combined difference between the simulation metallic- 
ities and the Zcmb estimate is a full order of magnitude; 
moreover, Figure [3] strongly suggests that the simulated 
metallicities are more likely underestimates than overes- 
timates, which would imply that an even larger fraction 
of stars formed with Z > Zqmb or Z > 10 Zcmb than 
has been calculated. However, even with this drastic 
change, half of the stars formed at 2.7 < z < 6 are formed 
in the CMB-regulated regime, a fraction that drops to 
10-15% at higher redshift. 

Our first major conclusion is therefore that if the CMB 
regulation hypothesis is correct, the CMB must have in- 
fiuenced the majority of star formation at 3 < z < 6, and 
may have been important to even higher redshift. 

4.2. Supernova and Gamma Ray Burst Rates 

Stars whose initial masses are at least Msn ^ 8 Mq 
end their lives as core-collapse supernovae. The fraction 
of stars above this limit, and therefore the number of 
supernovae produced per unit mass of stars depends on 
the IMF: a top-heavy IMF produces more supernovae. 
For a truncated power-law IMF, this ratio is given by 

iVsN ^ (a + 2)(AW"+^ ~MsN°+^) , . 

M* (a-f l)(M„iax°+2-Afmin"+2) ^' 

(|Bailin fc Harris! I2OO90 . The IMF may be top-heavy in 
the sense of having a flatter slope a, or be "bottom-light" 
in the sense of having a higher minimum stellar mass 
Mmin- To give an example of the quantitative size of the 
effect, flattening the slope from the canonical Salpeter 
value of a = —2.35 to a = —2.1 increases the supernova 
rate by 50%, while increasing Mmin from 0.30 Mq to 
5 Mq increases the supernova rate by a factor of 3.5. If 
the CMB-regulation hypothesis is correct, and the ma- 
jority of star formation at 3 < z < 6 occurred with a 
top-heavy IMF, t hen predictions of the hig h-redshift su- 
pernova rate (e.g. lDahlen fc FranssonlflQQQi ) are underes- 
timates by a corresponding factor. Some indirect conse- 
quences of this are discussed in sections 14.31 and 14.51 

The link between core-collapse supernovae and long 
duration gam ma ray bu rsts (GRBs) is now well estab- 
lished { Wooslev &: Bloom 2006). One might therefore 
expect that a top-heavy IMF, which increases the super- 
nova rate, might also increase the GRB production rate. 
However, there is both observational and theoretical ev- 
idence that GR Bs are pr eferentially produced by metal- 
poor stars (e.g. lYoon et al. 2006; Wolf & Podsiadlowskl 
120071 ). As CMB regulation resuhs in a top-heavy IMF 
only at high metallicities, it is not obvious whether it 
would result in an enhancement of the GRB rate with- 
out more detailed modelling. 

4.3. Reconstructing the cosmic star formation history 

A key goal in the study of galaxy formation is to recon- 
struct the cosmic star formation history (e.g. iLillv et all 
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ll996HMadau et al.lll996l ). There are three types of meth- 
ods for performing this measurement: the first is to di- 
rectly measure indicators of current star formation (such 
as Ha, UV, far-infrared, or radio continuum emission) 
in galaxies at a variety of redshifts; the second is to de- 
construct the stellar populations of low-redshift galaxies 
to determine their ages; and the third is to examine the 
stellar mass density in galaxies at a variety of redshifts, 
which is the integral of the star formation rate minus the 
stellar death rate over time. The first of these methods 
is direct: it measures star formation as it happens, while 
the latter two methods are indirect and rely on the stars 
that are produced. 

A key assumption in the indirect methods is the IMF: 
the mass of stars remaining after a given length of time 
from a fixed burst of star formation depends strongly on 
their mass distribution. If the IMF is very top-heavy, 
then the number of extant stars a given length of time 
later will be smaller than if the IMF is normal. In the 
extreme case where all stars are massive, then a burst of 
star formation will leave no visible stellar population a 
short time after the burst, and will be completely missed 
by any method that relies on the stars to measure the 
rate of star formation. 

Some of the direct star formation rate indicators also 
depend on the IMF. For example, most radio continuum 
emission in star-forming galaxies is synchrotron radiation 
from relativistic electrons that are accelerated in super- 
nova remnants (jCondonI Il992| ) . As discussed in § 14.21 
a top-heavy IMF increases the supernova rate with re- 
spect to the star formation rate, and therefore must also 
increase the radio continuum emission. Therefore, a top- 
heavy IMF both increases the derived star formation rate 
from direct indicators, and decreases the derived rate 
from indirect indicators. 

In fact, these measurement methods do not agree. At 
high redshift, the direct measurements of the star for- 
mation rate are systematical ly higher than the in direct 
methods. This can be seen in lHeavens et al.l ()2004[ ). who 
perform spectral synthesis modelling of SDSS galaxies 
at low redshift to determine their stellar populations, 
and use them to infer the star formation rate at red- 
shifts up to 10. They find that their star formation rates 
are systematically lower than the direct measurements in 
their higher redshift bins, precisely where CMB regula- 
tion would be expected to be impor tant. 

A similar result comes from Wilkins et al.l (|2008f ). who 
used measurements of the stellar mass density as a func- 
tion of redshift. The time derivative of this is the net 
change in stellar mass, which is the rate of star formation 
minus the rate that stars die. The latter term is a strong 
function of the IMF, and is much higher for a top-heavy 
IMF, requiring a much larger star formation rate to fit 
the same data. These authors find that their implied star 
formation history is consistent with direct measurements 
out to z = 0.7, but becomes increasingly discrepant at 
higher redshift, reaching a factor of 4 by z = 3. They 
suggest that an increasingly top-heavy IMF could ex- 
plain this discrepancy, but do not offer a physical mech- 
anism for this evolution. Although CMB regulation of- 
fers a natural explanation for why the IMF would ch ange 
with time, one problem is that lWilkins et al.l (|2008!D find 
evidence for a discrepancy between the direct and indi- 
rect star formation rates down to redshifts of z = 0.7, 



while CMB regulation is unlikely to operate much below 
z = 2.7. However, the discrepancy is relatively mild be- 
tween z = 0.7 and z = 2, and it is only beyond z = 2 
that the differences become significa nt. 

Similar considerations drove IDavel (|2008| ) to suggest an 
average IMF that slowly evolves with time; although he 
mentions the infiuence of the CMB as a potential driver 
for this evolution, the model is expressed as an overall 
gradual change with redshift rather than a change that 
specifically operates in the high-metallicity high-redshift 
regime. Fardal et al. (2007) perform a similar exercise 
while also considering the energy contained in the extra- 
galactic background light (which is dominated by light 
from stars formed at z < 2.6) as a constraint, and find 
that the best way to reconcile the measurements is with 
a top-heavy IMF, which they suggest occurs in the star- 
bursting galaxies that are increasingly common at red- 
shifts beyond where the extragalactic background is gen- 
erated. A CMB-regulated top-heavy IMF would perform 
a very similar function. 

4.4. Globular Cluster Bimodality 

A conundrum in the study of globular cluster systems 
around galaxies is the o bserved bimodality of their metal- 
licity distribution fe.g. iHarris et all 120061) . A possible 
explanation is that metal-poor clusters formed during 
the initial collapse of the protogalaxy, while the metal- 
rich population formed during subsequent major merg- 
ers. However, the metal-poor peak covers a narrow range 
in metallicity, which appears to be impossible to repro- 
duce unless the for mation of metal-po or clusters was 
abruptly truncated (jBeaslev et al.ll2002l ). 

However, if CMB-regulated star formation results in 
mostly or exclusively massive stars, there is another pos- 
sibility: the first round of globular cluster formation did 
occur over a wide range of metallicities from log Z/Zq ~ 
—2.5 upward, but the ones with logZ/Z© > —1 would 
have formed with a top-heavy IMF. The massive stars in 
the cluster would have died quickly, shedding the major- 
ity of their mass in the form of stellar winds and super- 
nova ejecta. If the massive stars made up a significant 
fraction of the cluster mass, as they would if the IMF 
were sufficiently top-heavy, this dramatic loss of bind- 
ing mass would result in the gravitational unbinding and 
dispersion of the cluster. 

In other words, the sharpness of the low- metallicity 
peak and the relative absence of intermediate-metallicity 
clusters near log Z/Zq — 1 may be entirely an artifact 
of a top-heavy IMF: higher metallicity clusters formed 
along with their lower- metallicity brethren, but did not 
survive to the present day. 

In the left panel of Figure [5l we have repeated the 
shaded and light shaded regions from Figure [4] that de- 
note the range in redshift and metallicity where a nor- 
mal IMF would be expected, while the unshaded regions 
are those where CMB regulation, if it occurs, would re- 
sult in a top-heavy IMF. We have overplotted the for- 
mation redshifts and metallicities of Milky Way globu- 
lar clu sters from the main sequence fitting performed by 
iMarfn- Franch ct al. ( 2009) . We have used the ages based 
on the I Potter et al.. (,2007 ) isochrones, and assumed that 
a normalized age of 1 corresp onds to an absolute age of 
12.8 Gyr, consistent with the iDotter et al.l ()2007D mod- 
els. These absolute ages have been converted to forma- 
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Figure 5. Metallicity versus f ormation age of Mil k y Wa y glob- 
ular clusters. Data come from iMarfn-Franch et aLl 1120091 ) . The 
relative ages have been converted into formation redshift by assum- 
ing (1) the WMAP 5 recommended cosmological parameter values 
((Hinshavir ct al. 2009), and (2) an absolute age scale of 12.8 Gyr 
and 12.1 Gyr for the left and right panels respectively. The dark 
shaded region is the region expected to exhibit a normal IMF, and 
is identical to the shaded region in Figure [T] The light shaded 
region also includes metallicities up to 1 dex larger, as in Figure|4] 



tion redshift by assuming the WMAP 5 recommended 
cosmological parameter values including the results of 
WMAP, baryon accoustic oscillations, a nd type la su- 
perno vae (the WMAP-hBAO+SN set of lHinshaw eFall 
l2009f ). While the relative ages are precise to within 2%- 
7%, at high redshift a very small variation in age corre- 
sponds to a large variation in redshift, and so the hori- 
zontal error bars in Figure [5| are very large at h igh red- 
shift. We have assumed the lZinn fc Westl ()1984[ ) metal- 
licity scale. 

Examination of the left-hand panel of Figure [5] reveals 
that there are many GCs that have metallicities signifi- 
cantly larger than the Zqmb at which they formed, i.e. 
above the dark shaded region, even given the large un- 
certainties in the formation redshifts due to the finite 
precision of the relative ages. However, if we assume 
that normal star formation can occur at metallicities up 
to 1 dex higher than ^cmb, i-e- in both the dark and 
light shaded regions, then the vast majority of GCs lie 
within the permitted region and only a few spill over into 
the forbidden region; the most discrepant clusters are 
Lynga 7, NGC 6171, NGC 6717, and NGC 6362. Are 
the existence of these clusters then conclusive evidence 
that CMB regulation is not important? 

It is interesting that three of these four clusters, 
NGC 6171, NGC 6362, and NGC 6717, lie in the less 
populated logZ/Z© ^ —1 metallicity regime between 
the metal-poor and metal-rich peaks, and that all of 
these clusters are quite low mass. Perhaps these objects 
are incompletely-disrupted remnants of a larger popula- 
tion of clusters that once inhabited the same parame- 
ter space, or perhaps the transition from a normal IMF 
to a strongly top-heavy IMF is gradual and these clus- 
ters formed with only a slightly top-heavy IMF. In these 
cases, we might expect to see systematic differences be- 
tween the stellar mass functions of these clusters com- 
pared to those of other Galactic GCs; detailed luminosity 
functions of these clusters would be well worth measuring 
to assess the viability of this explanation. 

However, we caution against overinterpretation of the 
location of these clusters. A serious concern is the ages: 
while the relative ages of clusters are precise, there is still 



large systematic uncertainty in the absolute age scale. 
The effects of this can be seen clearly by comparing 
the left-hand panel of Figure [5] to the right-hand panel, 
where we have reduced the absolute age scaling by 5% to 
12.1 Gyr. With this small change, the majority of GCs 
are consistent with forming at metallicities below Zqmb , 
and the only cluster that is clearly in the forbidden region 
is Lynga 7. The age of this particular cluster is suspect, 
however; other determinations of its age make it much 
younger (|Ortolani et al.lll993t ISarajedinil 120041 ) . and its 
location in the bulge introduces large uncertainties due 
to significant reddening. 

Because of the large inherent formation redshift un- 
certainties, we are not able to come to any firm conclu- 
sions regarding whether CMB-regulated star formation 
has shaped the population of extant globular clusters, 
although they appear to be consistent with the CMB- 
regulation hypothesis if the few discrepant clusters are 
indeed remnants of a disrupted population or simply ob- 
servational spill over. Perhaps if more precise absolute 
ages of GCs can be determined in the future, or if the 
detailed luminosity functions of the interesting clusters 
can be measured, it may be possible to address this issue 
more conclusively. 

4.5. Abundances and Abundance Patterns 

A top-heavy IMF in metal-rich high-metallicity star 
formation may leave characteristic patterns in the evo- 
lution of both the global metallicity and the abundances 
of particular elements. Not only does a top-heavy IMF 
increase the number of supernovae produced, but it also 
increases the average mass of the supernova progenitor^. 
More massive stars eject a larger fraction of their mass 
in the form of newly-synthesized metals, and are particu- 
larly efhcient at forming a elements such as oxygen, neon 
and magnesium (e.g. Wooslev fc W eaver 1995i). 

In fact, there is a slight inconsistency in using our simu- 
lations, which assume a normal IMF, to predict the frac- 
tion of star formation that occurs with a top-heavy IMF, 
as in §[3] If a large fraction of high-redshift star formation 
occurred in the CMB-regulated regime, as we predict, 
then they must have ejected a larger mass of metals into 
their environment, and therefore the global metallicity 
of the universe must have risen more steeply when CMB 
regulation was important, between at least z = 6 and 
z = 3, than the simulations predict. This would result 
in an even larger fraction of high-redshift star formation 
occuring in the CMB-regulated regime; however, we al- 
ready predict that the majority of high-redshift star for- 
mation occured above Zcmb, so this inconsistency does 
not qualitatively affect our conclusions. 

However, this more rapid rise in the global metallic- 
ity is precisely what is required to reconcile the discrep- 
ancy between the metallicities predicted by the simula- 
tions, log Z/Zq ^ —0.5 at z = 3, and those observed, 
which reach solar metallicity by the same redshift (see 
Figure [3]). 

A top-heavy IMF may also leave an imprint on the 
abundance patterns. Core collapse supernovae from 

^ At least, for an IMF that is truly top-heavy in the sense of 
having a flatter high-mass slope. This is not true of a "bottom- 
light" IMF, which is simply missing some fraction of low-mass stars 
but has otherwise the same functional form. 
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high mass progenitors eject a larger fraction of their 
mass in the forr n of g-elements, especiany ^^O, but also 
^°Ne ^ and ^-^Mg (jWooslev fc WeaveiifTOOalNomoto eFal] 
119971 ) . Therefore, a top-heavy IMF may produce more a 
elements than a population with a normal IMF. 

Current Galactic chemical evolution models accurately 
reproduce the observed abundance patterns of Galac- 
tic stars over a large range of metallicities, leaving lit- 
tle room for significant changes. However, one element 
whose abundances are not well reproduced is Mg, which 
is systematically underpredicted by the models over the 
range —2.5 < [Fe/H] < —0.5 (and slightly overpredicte d 
at lower metallicities; see figure 1 of lFrancois etal|[200l . 
Interestingly, these are precisely the metallicities where 
most of the CMB-regulated star formation occurs within 
our simulations, and ^^Mg would indeed be produced 
at greater rates with a top-heavy IMF. However, the O 
abundances do not support this conjecture: they are, if 
anything overpredicted at these metallicities, while a top- 
heavy IMF would have an even larger effect on ^^O than 
on ^^Mg. Unless some unaccounted-for process system- 
atically drives down the ^^O abundance by a significant 
amount, a more parsimonious explanati on is that the 
yields require adjustment, as suggested bv lFrancois et al.] 
((200411 . 

5. CONCLUSIONS 

We have inves t igated the suggestion of ILarsoiil (|2005[ ) 
and lSmith et al.l ()2009[ ) that star formation at high red- 
shift and high metallicity results in a top-heavy IMF due 
to the infiuence of the CMB as an effective tempera- 
ture floor for the collapsing protostellar core. By assum- 
ing that CMB regulation is important if the minimum 
temperature the core would reach in the absenc e of th e 
CMB, according to the models of lOmukai et al.l ()2005[ ). 
is less than the CMB temperature at a given redshift, 
we are able to parametrize the evolution of the critical 
metallicity Zqmb as a function of redshift. Comparison 
with t he 3D hydrodynamical simulations of ISmith et al.l 
(|2009D suggests that our model is generally accurate, but 
may underestimate Zqmb by up to an order of magnitude 
in particular cases. 

By comparing Zqmb to the metallicities of stars formed 
in high resolution cosmological simulations of the forma- 
tion of eight ^ galaxies, we conclude that if CMB reg- 
ulation does operate, the majority of stars at redshifts be- 
tween 2.7 and 6 would have formed in the CMB-regulated 
regime. We have investigated five possible observational 
signatures of CMB regulation: 

1. Current predictions of supernova rates, and possi- 
bly GRB rates, are underestimates at z > 3. 

2. Indirect measures of the star formation history 
based on the stellar populations left behind by 
episodes of star formation, such as the age distri- 
bution of local stellar populations and evolution of 
the stellar mass density as a function of redshift, 
would systematically underestimate the star for- 
mation rate at redshifts greater than 2.7 compared 
to measurements based on UV and far-infrared fiux 
generated directly by star formation. Such a dis- 
crepancy is seen, although it appears to extend to 
lower redshift. 



3. Globular clusters that formed at high metallicity 
and high redshift may have quickly evaporated, 
leaving no present-day evidence of their existence, 
possibly explaining the relatively narrow metallic- 
ity range of the metal-poor GC population. In this 
case, we would expect there to be no extant GCs 
that are both metal-rich and formed at very high 
redshift. Because of the difficulties in determining 
absolute ages to old GCs, and the steepness of the 
age-redshift relation at high redshift, current data 
are not able to constrain this hypothesis. 

4. Cosmic metallicity may rise more rapidly at early 
times than predicted by current simulations that 
assume a normal IMF, reaching solar metallicities 
at higher redshift. This could reconcile the simula- 
tions with observations that show a mean metallic- 
ity of approximately solar for stars formed out to 
redshift 3. 

5. Elements produced by high-mass supernova pro- 
genitors, particularly the a elements O, Ne, and 
Mg, would be produced with greater abundance 
between redshifts 6 and 3, or [Fe/H] ~ —2 to —0.5, 
than predicted in models that assume a normal 
IMF. This could help explain the high Mg abun- 
dances seen in stars of this metallicity range in the 
Galaxy, but would exacerbate the overprediction 
of O over the same range in metallicities. Since 
the O yields are more sensitive to the IMF, we 
conclude that the abundance patterns are not pro- 
duced by a top-heavy IMF, although a bottom-light 
IMF would not alter the abundances. 

In conclusion, the observations neither provide conclu- 
sive evidence for the hypothesis that the CMB can reg- 
ulate star formation, nor rule it out. The most convinc- 
ing evidence would come in the form of either the clear 
presence or absence of an envelope in redshift-metallicity 
space, above which stars are either absent or have greatly 
reduced numbers compared to the number of stars be- 
lieved to have formed there from other measurements; 
unfortunately, as discussed in § 14.41 it is unclear if stel- 
lar ages will ever be known precisely enough to perform 
this test. Another potential avenue of exploration would 
be to use the methods that currently argu e for a mean 
cha nge in the IMF with time, such as in [Pavl ([2008h 
and iWilkins et all ()2008[ ). but subdividing the popula- 
tions by metallicity to determine if an evolving IMF is 
only required for metal-rich stars, as would be expected 
if the evolution in the IMF is due to CMB regulation, or 
if it must be truly universal. 

This paper makes use of simulations performed as 
part of the SHARCNET Dedicated Resource project: 
"MUGS: The McMaster Unbiased Galaxy Simulations 
Project" (DR316, DR401, and DR437). 
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